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COMBUSTION OF LIQUID SPRAYS AT HIGH PRESSURES 


Summary 

This report discusses activities under NASA Contract NGR 39-009-077 
for the period January 1, 1976, to December 31, 1976. During this period 
the combustion of fuel sprays in a stagnant environment was considered 
both theoretically and experimentally. The work included operating 
conditions where liquid fuel approached its thermodynamic critical point 
during combustion. 

The experiments considered methanol and n-pentane as fuels, injected 
through a single-hole, orifice-type injector into pure air. Test pressures 
were in the range 0.1-9 MPa. Measurements were made of spray and flame 
boundaries. The experiments were compared with theoretical predictions 
based on a locally homogeneous two-phase flow model. The turbulence charac- 
teristics of the jet were represented by an integral model, using a variable 
density entrainment law which had been developed for gas-gas and gas-liquid 
jet processes. The theory had not previously been compared with sprays. 
Aside from spray and flame boundaries, the model estimates profiles of 
mean quantities within the spray. 

The theory correctly predicted the trends of the data, but generally 
underestimated the extent of the spray and flame boundaries by 30 - 50 per- 
cent. The results indicate that slip effects were still important for 
the present experiments (the Sauter mean diameters of the sprays were 
approximately 30 ym at atmospheric pressure under cold flow conditions) . 

The accuracy of the predictions is poorer at high pressures, even though 
the density ratio of the two phases approaches unity, which should improve 
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the locally homogeneous flow approximation. The sprays are shorter at 
high pressures, however, and significant slip effects near the Injector 
appear to override the density ratio effect. 

The predictions indicate regions within the spray where water vapor 
produced by combustion should condense, however, the condensed water 
boundary always fell within the spray boundary for the present test 
conditions. Both theory and experiment did not indicate any unusual 
phenomena when spray gasification was completed by the fuel passing 
through its thermodynamic critical point. 

The locally homogeneous model developed in this study is convenient 
to use and requires a minimum amount of input data. Quantitative accuracy 
could be improved by adjusting empirical parameters in the model from 
the values that were optimized for gas-gas and gas-liquid jets, however, 
further data over a wider range of injector conditions would be desirable 
prior to such recorrelation. 
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ABSTRACT 

The combustion of a liquid fuel spray produced By a plain 
orifice pressure atomized injector was considered. Combustion was 
examined in a stagnant air environment at pressures of .1—9 MPa for 
n-pentane and methanol. The higher pressure levels are in excess 
of the critical pressure of the fuels. Measurements of spray and 
flame boundaries were compared with, predictions based on a locally 
homogeneous model of the flow which had been developed for turbulent 
gas-gas and gas-liquid j ets . The theory correctly predicted the 
trends of the data, but underestimated the extent of the spray and 
flame boundaries by 30-50 percent. The results indicate that slip 
is important for the sprays considered in this investigation (cold 
flow Sauter mean diameters were approximately 30 ]im) . No unusual 
effects were observed for supercritical conditions. A region where 
water vapor produced by combustion-should condense was found 
theoretically, however, the water condensation region fell within 
the spray boundaries for the present test conditions. 



CHAPTER I 


INTRODUCTION 

1.1 General Statement of the Problem 

Many Important combustion systems such as liquid propellant 
rocket engines, diesel engines and gas turbines depend upon the spray 
combustion process. Modem design of combustion devices has emphasized 
liquid fuel combustion at high pressures often approaching the critical 
point of many fuels. Therefore an understanding of the spray 
combustion at near critical conditions would provide a useful design 
tool for combustion chamber development. 

The spray combustion process consists of a liquid fuel flowing 
into a gaseous enviroiunent and reacting with the surrounding gas . 

Figure 1 illustrates the regions present in a spray combustion process 
as well as typical velocity and concentration profiles. All of the 
injected liquid is contained within the spray boundary. Within this 
region the liquid jet leaving the Injector breaks up into droplets. 

The droplets are heated by the gas in the region and evaporate 
supplying gaseous fuel. The fuel is then transported toward the 
reaction zone by turbulent mixing. The reaction zone is a region 
where the oxygen and the fuel are present and the chemical reaction 
occurs. The reaction zone is defined by the position where the 
mean oxygen concentration goes to zero on the inside and where the 
mean fuel concentration vanishes on the outside. Combustion products 
produced in the reaction zone are mixed throughout the flow. 
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In order for the jet to spread, additional fluid is added to the 
flow by an entrainment process. Oxygen and all nonreactive gases are 
added to the flow in the same ratio as they present in the surrounding 
fluid. The jet boundary identifies the outer region of the flow where 
the velocity goes to zero. 

The typical variation of the centerline temperature is shown 
in Figure 2. The centerline temperature remains relatively low in 
the region where drops are present. The droplets gasify by a 
conventional evaporation process at low pressures; at high pressures 
gasification can occur by the drops exceeding their thermodynamic 
critical point. The centerline temperature increases from the spray 
tip toward the reaction zone where a maximum value is reached. As 
additional fluid is drawn into the flow, the temperature decreases 
downstream of the flame tip, similar to processes occurring in a 
nonreactive jet. 

The present study concentrated on examining the spray combustion 
process at elevated pressure conditions typical of spray combustion 
devices. The spray boundaries and the flame shape were determined for 
several fuels at near critical and supercritical conditions. 

1.2 Previous Related Studies 

The combustion of a liquid fuel spray is a complex phenomenon 
requiring knowledge of the evaporation and combustion processes. 
Individual droplet combustion has already been considered in some 
detail (1~4) . These studies have provided the basis for predicting 
gasification rates typical of spray combustion systems. 
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In order to utilize the individual droplet studies, the structure 
of the spray must be identified to specify the local environment of 
the droplets. Only recently have Investigations of the spray structure 
during combustion appeared due to the experimental difficulty of this 
type of study. However, the study of some nonreacting sprays provides 
some insight into the two-phase nature of the spray combustion process. 

Hetsroni and Sokolov (5) investigated the effect of very small 
liquid droplets on the structure of turbulent air jets. The droplet 
distribution and velocity profiles were measured and correlated with 
a Gaussian distribution. It was concluded that the velocity profiles 
of single-phase and two— phase jets were similar. 

Weimer, Faeth and Olson (6) demonstrated that two-phase 
turbulent jets of gases into liquids may be modeled together 
irrespective of the vapor liquid system. Vapor penetration lengths 
for condensing water, ethylene glycol and iso-octane jets were 
correlated using a variable density single-fluid model for the 
two-phase flow and a variable density entrainment law to account 
for the turbulent mixing of the injected fluid and the surroundings. 
Results from Kemey, et al., (7), for condensing steam-water jets, 
were also correlated by the variable density model. 

Tross (8) experimentally investigated a turbulent two-phase 
air-water jet. It was concluded that the velocity, void fraction 
and momentum flux possessed similar profile distributions at any 
position in the flow when the ratio of these quantities to their 
centerline values was plotted against a nondimens ional radial 
coordinate. Gaussian distributions were shown to fit these profiles. 
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The characteristics of nonreacting diesel fuel sprays were 
experimentally investigated by Wafcuri, et al. , (9). Photographic 
measurements of fuel sprays were recorded using a high-speed motion 
picture camera. Nondimensional correlations were developed for 
spray length as functions of density ratio, inlet velocity, initial 
diameter and time. It was concluded that for well-atomized sprays 
there was negligible relative velocity between the droplets and the 
surrounding fluid and the two phases could be treated as a homogeneous 
mixture. 

Newman and Brzustowski (10) investigated the behavior of a 
turbulent two-phase jet near the critical region of the injection 
fluid. Liquid carbon dioxide was injected into stagnant gaseous 
atmospheres of carbon dioxide and nitrogen. The two-phase spray was 
treated as a homogeneous turbulent jet by assuming that the droplets 
moved at the same velocity as the surrounding gas and were in thermal 
equilibrium with the gas. Success of the model was limited to order 
of magnitude predictions on the size of the spray boundary. 

Similarities also exist between two-phase reacting sprays and 
single-phase reacting jets. Thrlng and Newby (11) analyzed the spray 
combustion process by burning atomized oil jets. The combustion 
length of the oil jets was primarily governed by the effectiveness 
of the mixing process. It was shown that the mixing process for two- 
phase jets could be adequately described by a mixing model developed 
for nonreacting flows. 

Further investigations of the spray combustion process have been 


performed under atmospheric pressure conditions by Chigier and 
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coworkers (12, 13) and Onuma and Ogasawara (14). In th,ese 
investigations, the spray structure was defined by measuring gas 
temperatures, droplet sizes j velocities and concentrations using air 
atomizing nozzles. The studies indicate that droplets in a flame do 
not bum individually but that the fuel vapor diffuses from the 
region of drop evaporation and bums like a gaseous diffusion flame. 

Khalil and Whitelaw (15) also investigated the two-phase spray 
combustion process, but found somewhat different results. 

Experimental values of velocity, turbulence intensity, temperature 
and droplet concentration were determined for a hollow cone kerosene 
spray at atmospheric pressure. The flame length increased 
significantly when the Sauter mean diameter of the spray was increased 
from 45 pm to simulated gaseous diffusion flame which predicted more 
rapid development of the flame. This suggests that the assumption 
of locally homogeneous flow, which does not provide for drop size 
effects, can be invalid even for sprays having a Sauter mean diameter 
as small as 45 pm. 

Avery and Faeth (16) investigated a much different problem 
involving the combustion of a gaseous oxidizer jet into a liquid metal 
fuel. The relative velocity between the two phases is small because 
of the relatively low inertia of the gas compared to the liquid fuel. 
Density variations were handled by a coordinate transformation 
reducing the system to the case of an incompressible jet. A unified 
correlation of flame length, temperatures, and velocities was 
developed for both the two-phase system and earlier studies of 
single-phase reacting jets. Figure 3 presents an example of the flame 
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Figure 3 Comparison of Predicted and Measured Penetration Lengths 
for Gas-Liquid and Gas-Gas Turbulent Jets, Reference (16) 



9 


length results. The gas-liquid system, of Reference (16) was compared 
with the gas-gas results of Hawthorne, et al., (17) and Wohl, et al., 
(18) . It was also possible to correlate the measurements of 
condensing vapor jets, as shown in Figure 3. 

Under the proper conditions the structure of two-phase turbulent 
jets are similar to single-phase jets and more importantly two-phase 
spray combustion processes are similar to gaseous diffusion flames. 
This is particularly the case when the droplets are small enough 
that slip between the two phases is negligible. 

1.3 Specific Statement of the Problem 

The preceding discussion has indicated that although some work 
has been done on the structure of the two-phase turbulent spray 

r 

combustion processes, few studies have examined the spray combustion 
process at elevated pressures. Under these conditions the density 
variation between the liquid and g aseou s phases becomes small, 
implying a closer approximation to the no-slip model. 

The present study will examine the injection of a liquid fuel 
spray from a single-hole orifice-type injector, without any swirl, 
issuing into a stagnant environment of pure air. With this in mind 
the objectives of the present study are as follows: 

1) Experimentally determine the spray boundary of a two-phase 
nonreacting liquid fuel jet. 

2) Experimentally determine the spray boundary and the flame 
shape of a reacting two-phase turbulent jet, 

3) Compare the predictions of the two-phase spray combustion 
model proposed by Avery and Faeth (16) to the experimental 
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results for reacting and nonreacting liquid fuel 
sprays. 



CHAPTER II 


EXPERIMENTAL APPARATUS AND PROCEDURE 

2.1 Description of the High Pressure Experimental Apparatus 

The general purpose of the experimental apparatus was to provide 
a large quiescent high-pressure environment for observation of the 
fuel spray- High-speed photographs of the spray were taken at several 
axial locations to form a complete picture of the spray. The 
photographs were used to determine the spray boundaries and flame 
shapes. A schematic representation of the overall experimental 
apparatus is illustrated in Figure 4 and a photograph of the facility 
is shown in Figure 5. 

3 

The test chamber consists of a 9000 cm cylindrical vessel 66 cm 
long with an internal diameter of 13 cm. Maximum working pressure of 
the chamber was limited to 10 MPa. —The vessel was constructed from 
low carbon steel which necessitated coating the interior surfaces with 
a rust- inhibiting paint to prevent corrosion by the combustion 
products. The test section of the chamber was located in the upper 
portion of the vessel in order to maintain the tip of the flame as 
far as possible from the end of the chamber. Photographic observation 
of the fuel spray was performed through two quartz windows 
approximately 1.25 cm in diameter located on opposite sides of the 
vessel. Visual observation of the spray was conducted through a third 
quartz window. 

Compressed air for the experimental apparatus was supplied by a 
reciprocating type Ingersol Rand Compressor capable of supplying air 
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up to 20 MPa. The compressed air was filtered at two locations by 
Matheson type 451 and 453 filters to remove oil, water, or other 
impurities to the 5 ym level. Chamber pressure was monitored by one 
of three Heise pressure gages depending on the operating pressure. 

Tlie three pressure ranges were 7500, 1500 and 300 psi with a specified 
accuracy of 7.5, 1.5 and .3 psi, respectively. 

A sketch of the fuel injection assembly is shown in Figure 6. 

The assembly consists of a fuel loop, nozzle, nichrome igniter coil 
and alignment supports. The loop configuration was selected to store 
the fuel in order to maintain a continuous fuel slug prior to 
injection. Fabrication of the loop used one-quarter inch stainless 
steel tubing. The nozzle used throughout the tests was a stainless 
steel straight hole type injector, model .000009 solid stream tip, 
supplied by Spraying Systems Company. Attachment of the nozzle to 
the fuel loop was accomplished by welding the nozzle to a one-quarter 
inch Swagelok union which was connected to the fuel loop. The 
position of the nozzle was checked prior to each test by examining 
the fuel Injection assembly through the camera eyesight. The 
alignment supports were also used to position the ignitor coil. 

Injection of fuel into the chamber was accomplished by 
applying a differential pressure across the fuel loop forcing the fuel 
slug through the nozzle. The driving pressure was adjusted by a gas 
regulator and monitored by a differential pressure gage comparing the 
driving pressure to the chamber pressure. A solenoid valve was 
located upstream of the fuel loop in order to rapidly apply the 
driving pressure for injection. A 500 ml capacitance chamber was 
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Figure 6 Sketch, of the Fuel Injection Assembly 
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located in series with, the solenoid valve in order to maintain the 
driving pressure relatively constant during the injection process. 

The test procedure was designed to use only a small amount of fuel 
in order to avoid large changes in the chamber pressure, temperature, 
and gas composition during a test. The total quantity of fuel injected 
was limited to 4 ml for all tests. 

Ignition of the fuel spray was accomplished by positioning a 
heated coil of nichrome wire near the fuel spray. The coil was 
approximately 2.5 cm long and 0.5 cm in diameter constructed from 28 
gage wire. The power input to the coil was regulated by a variable 
transformer and increased until the coil began to glow red. 

Examination of dark field photographs of the near injector region 
confirmed that the flame was attached back to the nozzle. 

Photographic measurements of the fuel spray were recorded with a 
Photosonic 1-B motion picture camera using a Kepco SM36-5 AM d.c. 
power supply. The film speed was indicated by a timing light on the 
camera activated by a Wollensak Pulse generator, model 3106A, set 
at 100 pulses/s. Backlighting for the shadowgraph spray measurements 
was supplied by a Pek, model 401A arc lamp using a 75 watt mercury 
bulb. The light was focused into a parallel beam using the optics 
located in the arc lamp. A diffuser screen was used to equalize 
the intensity of the light beam. Kodak plus-X reversal film was used 
for all photographs. 

The entire injection process was controlled by a mechanical 
timer. At the beginning of the timing cycle the electrically driven 
motion picture camera and pulse generator were activated and the 
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camera was permitted to reach, operating speed. The solenoid valve 
was then actuated driving the fuel into the chamBer. Injection 
continued for approximately 3 seconds before the solenoid valve was 
closed. The camera power and timing light generator were then turned 
off to complete the test cycle. 

2.2 Description of the Low Pressure Experimental Apparatus 

The purpose of the low pressure apparatus was to measure the 
flame boundaries at atmospheric conditions. A sketch of the apparatus 
is shown in Figure 7. The fuel injection and ignition systems were 
identical to those used in the high pressure tests with the exception 
that the nozzle was oriented horizontally. Photographs of the spray 
were recorded with a Graphlex camera. Polaroid type 57 film was used 
in these tests. 

2.3 Description of the Experimental Procedure 

The present investigation considered both methanol and n-pentane 
fuel sprays at ambient pressures of .1, 3, 6 and 9 MPa. The fuel 
used was of certified grade supplied by Fisher Scientific Company 
with a minimum purity of 98 percent for the pentane and 99.9 percent 
for the methanol. The temperature of the high pressure environment 
and liquid fuel was taken to be at room temperature which varied from 
20°C to 30°G. Spray boundaries were recorded for both the reacting 
and nonreacting sprays along with the flame boundaries for the 
reacting spray. 

2 

The field of view consisted of an area approximately 1.25 cm . 
Because the fuel spray xras longer than the field of view, the Injector 
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Figure 7 Sketch of the Low Pressure Apparatus 
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had to be moved In order to obtain a complete picture of the process. 
For each location several tests were recorded in order to obtain an 
average value for the boundaries at a particular test condition. 

At each test pressure, shadowgraph photography was used to 
determine the spray boundaries; dark field photographs were ^ployed 
for the flame boundaries. A typical example of a turbulent spray 
flame is shown in Figure 8. The fluctuation in the boundaries due 
to the unmixedness of the spray process are very apparent. The 
irregularities in the boundaries demonstrate the need for a 
statistical approach in determining the boundaries. 

The photographs were analysed frame-by-frame using a microfilm 
viewer. The position of the boundaries was determined by measuring 
the distances on the viewer and then scaling the measurement back 
to true size based on the nozzle width. Radial positions were measured 
at several axial locations along the spray and the spray penetration , 
distances were also recorded. These distances were then averaged in 
order to determine the mean positions. 

Each test required the pressure vessel to be disassembled. During 
this period the chamber was purged of any combustion products from 
previous tests and the windows cleaned of any soot and condensed water 
droplets. Fuel could then be loaded into the fuel loop and the 
chamber reassambled for the next test. 

The chamber was then pressurized to the given test condition. 

The pressurization process occurred very slowly in order to permit 
the gas and the fuel slug to come into thermal equilibrium with the • 
chamber. Once the given test condition was reached, the variable 




Figure 8 Photograph of a Pentane Spray Flame at Atmospheric Pressure 
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transformer was adjusted until tlie ignitor coil just Began to glow 
red and tB.e cycling timer actuated. 



CHAPTER III 


THEORETICAL CONSIDERATIONS 


3.1 Introduction 

The theoretical objective of this investigation was to examine 
the ability of an integral model, which assumes a locally homogeneous 
flow, to predict spray boundaries and flame shapes at high pressures. 
Major emphasis was placed upon determining the effectiveness of the 
model as the ambient pressure increased. The present theory is similar 
in many respects to that of Reference (16) for the combustion of a 
submerged oxidizer jet in a liquid metal. The major difference is 
that the present investigation must consider the evaporation of the 
liquid fuel droplets and the gas in the continuous phase. 

The general approach of the analysis will involve the assumption 
of similar profiles for various quantities. Integration of the 
conservation equations, in conjunction with an entrainment expression 
proposed by Morton Cl9)j yields correlations of enthalpy decrement, 
concentration, and velocity throughout the spray. The analysis also 
provides an estimation of the pentration length and the radial 
boundaries of the spray as well as the flame boundaries. 

The present investigation considers a reacting spray combustion 
process. The flow consists of a region of liquid fuel droplets 
surrounded by a region of fuel vapor and product gas, finally 
bounded by a reacting Interface with a stagnant air environment. The 
major assumptions in the analysis are as follows: 
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1) The flow is considered to be a steady axisymmetric Jet 
issuing into an infinitely large quiescent environment. 

2) The flow field is considered to be at constant pressure. 

3) The vapor and liquid phases are in local homogeneous 
equilibrium, i.e., they have the same local velocity and 
temperature. 

4) Only turbulent transport processes are considered. 

5) Radial gradient terms are much greater than axial 
gradient terms. 

6) Magnitudes of fluctuating quantities are small compared 
to magnitudes of mean values . 

7) The jet entrains material from the environment according 
to the entrainment expression developed by Morton C19) . 

8) Combustion occurs as a one-step chemical reaction. 

9) The reaction rate is infinitely great and the reaction 
. itself is localized in the flame zone. The effect of 

unmixedness is handled using the same approach as in 
Reference (16) . 

10) Buoyancy forces are neglected, 

11) Flow properties are Toniform across the width of the jet at 
the injector exit. 

12) Mean profiles of axial mass flux,, momentum flux and energy 
flux are assumed to be similar in shape at all positions 
in the jet. 

13) The ideal gas law will be used to model the gas phase 
properties. 
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14) Constant specific Bleats will Be employed to describe ttie 
enthalpy variations. 


3.2 Prediction of tli,e Penetration Length 

Under these assumptions the integral equations of conservation 
of mass, momentum, energy, and species are: 


dz 

OO 

X 

purdr = - (rpv)^^^ 

(3.1) 

2tt 

00 

[ 

2. 2 2* 
pu rdr = irr^ p^u„ = M 

^ O ^0 0 0 

(3.2) 

2tt 

Jo 

00 

f 

puAhrdr = m Ah 

(3.3) 

2ir 

Jo 

00 

f 

0 O 

puAyrdr = m Ay 
0 o 

• (3.4) 

where 

J” 



Ah 

II 

I 

(3.5) 


and 


Ay 


c y 

00 


yox 

+ 1 

y 

0X„ 


(3.6) 


The variable Ay defined in Equation (3.6) is a Shvab-Zeldovich 
variable which results from eliminating the reaction rate terras 
between the equation of conservation of species for the fuel and 
oxidizer. 

The initial condition on the above set of equations is 


CO 

2ir J purdr = m^ , z = 0 (3.7) 

The boundary conditions for a quiescent environment having a 


constant composition are: 
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u = Ah = Ay=0, P^Poo* r - (3.8) 

The dependent variables in Equations (3.1)-(3.4) are now 
normalized in terms of their centerline values at each axial station 
of the jet; 


f , (n) = pu/p u 
J- c c 

(3.9) 

f2(n) = pu^/p^u^2 

(3.10) 

f3(n) = puAh/Pj,u^Ah^ 

(3.11) 

f^Cn) = puAy/p^u^Ay^ 

k 

(3.12) 

r\ = r/a(z) 

(3.13) 


It is now assumed that the f^ are similar at all axial locations 
and the following integral constants are defined 

CO 

I. = r f .CTi)ridn ; 1 = 1 , 2 , 3, 4 (3.14) 

Jo 

The values of these integrals only change when the assumed profiles 
are changed, and they are constant for fully developed flow under 
the present assumptions. 

The entrainment of the jet is represented by a model suggested 
by Morton (19) , which has been found to be valid for the self- 
preserving region of the jet. The form of the entrainment expression 
is : 

1/2 

C- 2TTrpv)^ = 27Tp^Eau^(p^/pJ (3.15) 

where E is the entrainment constant for the jet. The value of the 
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entrainnient constant is dependent upon the assttmed shape of the 
profiles, the characteristic radial scale factor and the normalizing 
values of density and velocity used in the expression. 

The form of the equations may now be simplified by introducing 
a new radial scale suggested by Morton (19) for large density 
differences. The new scale length is defined by: 

b(z) = a(z) (3.16) 

and reduces the governing equations to a form that is similar to a 
constant density homogeneous jet. 

Substituting Equations (3. 9) -(3. 16) into Equations (3. I)- (3. 14) 
and rearranging yields the following set of equations. 


~ (u b^I ) = 
dz c i . 

Ebu^ 

(3.17) 


M 

o 

(3.18) 

2irp^,j_.b2Ah^l3 

= m h 
o o 

(3.19) 

2 

2irp u b Ay I, 
c -^c 4 

= m y 
o o 

(3.20) 


Transforming the initial condition yields 


2irp u b^I, = m 
“cl 0 


z = 0 


The solution of Equations (3, 17) -(3. 20) yields 


(Trp„M^) 


m 


o 


M, 



c 




(3.21) 


( 3 . 22 ) 


(3.23) 
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27Tl3^p«u^b2 


(3.24) 


(3.25) 


(3.26) 




where 


# = 1 + E 


r -] 1/2 (irp M )^^^z 
2 1 CO o 

^2 


(3.27) 




In order to apply these results to estimation of flame 
penetration length, the flame position must be specified. In the 
case of homogeneous laminar diffusion flames, the flame zone is very 
thin and it is reasonable to assume that the concentrations of fuel 
and oxidizer are zero at the flame. This provides a criteria Ay=l for 
locating the position of the flame. 

For turbulent flow, the flame zone is relatively thick and 

profiles of mean fuel and oxidizer concentrations overlap to a 

considerable extent (11, 17). Therefore, the outer limit of the 

combustion zone, corresponding to the point of disappearance of the 

injected material, occurs at values of Ay that are less than unity. 

In order to allow for this effect, it is assumed that the injected 

material is absent in the region defined by Ay ^ where is 

an empirical parameter. With this specification the maximum length 

of the combustion process, L , is determined by the criteria 

r 


Aye = 5 


z = L 


(3.28) 



28 


Applying the criteria of Equation (3.28) to Equations (3.25) and 
(3.26) yields the following expression for the penetration length 



(3.29) 


(3.30) 


A penetration length correlation for an evaporating jet can also 
be obtained from Equations (3.24) and (3.26). For an evaporating jet, 
the criteria for the disappearance of liquid is that the local 
enthalpy must be greater than the saturated vapor enthalpy at the 
chamber pressure. Allowing for the unmixedness of turbulent flow, 
as before, provides the following specification for the maximum 
lengthwise position of the liquid. 


Ah = £ (h - h ) = e Ah (3.31) 

C g g 00 g g 

where h is the saturated vapor enthalpy and e is an empirical 
8 g 

parameter less than unity. 

The correlation for liquid penetration length becomes 



r I, 1 

ri-1 

1 

Ll 

e Elo 

[2 J 

L g 3J 



1/2 


(iTpJ[^) 



d B 
o g 


where 


B 

g 



(3.32) 


(3.33) 
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The combined correlation of spray and flame lengths from Reference 
(16) is shown in Figure 3. 

A correlation of velocity along the centerline of the jet may 
also be obtained for large z/d^, by considering Equations (3.23) and 
(3.27). The centerline velocity may be represented by 


M 


* u 
o c 


^2 f2 ']1/2 -E(ttp«M^)^^^z 


^1 1 ^ 2 . 


(3.34) 


m^ 


This correlation is shown in Figure 9 from Reference (16). 

Avery and Faeth (16) demonstrated that enthalpy decrement ratios 
and nozzle fluid concentrations varied in the same manner, these 
results are illustrated in Figure 10. This suggests that and 
are equal. Equations (3.24) and (3.25) may be combined into the 
form 


Ah Ay 

c _ •'c 

The enthalpy decrement ratio may be calculated from the 
expression 



n 


m 


1 - i 

R m 

1 I 


(3.35) 


(3.36) 


where 


AT = T - T^ (3.37) 

The fluid concentrations may be derived from the Shvab- 


Zeldovich variables. 
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Figure 9 Correlation of Centerline Velocity, Reference (16) 
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Figure 10 Correlation of Centerline Scalar Quantities, Reference 
(16) 
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(3.38) 

(3.39) 

(3.40) 

(3.41) 


The temperature at the penetration length of the spray 
boundaries may be calculated using Equations (3.36)-(3.41) . 
Concentrations of the species are calculated at selected values of 
Ay/Ay and used to calculate Ah /Ah until the selected value of 

O CO 

Ay/Ay^ equals the calculated value of Ah/Ah^ as required by Equation 
(3.35). 


The penetration length for the fuel spray boundaries is 

calculated using the above procedure. The driving potential for 

evaporation B , in Equation (3.33) may be calculated using the 
S 

computed value of Ah/Ah^. The fuel spray penetration length in 
Equation (3.32) can then be determined. 

Analysis of the temperatures in the fuel spray region indicated 
that condensed water vapor is present in some circumstances. The 
calculation procedure can also be used to determine the spray 
penetration length of the condensed water vapor using Equations. 
(3.32) and (3.33). 

The presence of both liquid and gaseous fuel must be considered 
in the spray region. The vapor pressure of the fuel was calculated 
from a relationship of the form 




= A - 


B 

(C+T) 


(3.42) 
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The mass fraction of fuel may Be computed from the following relation 
i ^ 1 ^i 


^^F,G 




(3.43) 


i ^ 1 1 i 

The mass fraction of the liquid fuel can also be determined 


^F, L = ^F “ ^F^G 


(3.44) 


The enthalpy decrement in Equation (3.33) can be calculated treating 
each phase of the fuel as a separate component. 


3.3 Prediction of Radial Boundaries Accounting for Variable Density 
The position of the radial boundaries may be predicted by 
extending the analysis derived for the penetration length. Equations 
(3.25) and (3.27) are combined to predict the axial position for a 
given concentration decrement into the form 


z 

d 


Ay 

1 ■'o 


- 1 


1 1/ ^y 
L 4 


2 E 


2j]l/2 

^2. 


P "1/2 


Lp 


(3.45) 


The characteristic length scale may be computed by combining Equations 
(3.16) and (3.22) into the following relationship. 


a = 


$(12/2) 


2 )^/^ p 1^/2 

J Poo. 


(3.46) 


where 


$ = 1 + 2 E 


LI2J 


1/2tp -tl/2 


*^oJ 


(3.47) 


Tross (8) demonstrated that the normalized functions described 


by Equations (3. 9) -(3. 12) may be written as 
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f^(n) = Exp[- i = 1, 2, 3, 4 


C3,48) 


Substitution of Equation (3.48) into Equation C3.14) yields 

CO 

J i 2, 3, 4 (3.49) 

Integration of the above expression yields 


I = _JL ; i = 1, 2, 3, 4 (3.50) 

^ 2K. • 

1 

The nondimens ional functions represented by Equation C3.48) can now 
be expressed by the following relationship 


fjCn) 



(3.51) 


The dimensionless radial distance ri can be represented as a function 
of the concentration decrement by combining Equations C3.9), C3.12) 
and C3 . 51) 



By combining Equations (3.13), (3.16), 
distance may be expressed as 


(3.52) 

(3,46) and (3.52) the radial 



(3.53) 


3.4 Prediction of Incompressible Radial Boundaries 

The outer location of the liquid fuel droplets is bounded by the 
spray boundary while the outermost location of fuel vapor is 
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indicated by the flame boundary. It is necessary to determine a 
relationship between the radial spread of the Jet and axial position. 
The spreading rate of the jet will be taken to be the same as an 
incompressible jet. 

Comparison of the normalized radial boundaries measured by 
Tross (8) , Hetsroni and Sokolov (5) and McCreath and Chigier (13) 
indicate that the radial spread of the velocity in the spray is not 
strongly influenced by the density ratio of the spray to the ambient 
environment. Figure 11 compares the radial spread of steam-water 
jets, isothermal air jets and isothermal heterogeneous jets. These 
data encompass jets having density ratios greater than, equal to, and 
less than the ambient fluid. 

The influence of chemical reaction on the radial spread was 
considered by Chigier and Roett (12) . The effect of chemical 
reaction on the radial spread of homogeneous and heterogeneous jets 
is shown in Figure 12. The reacting and nonreacting heterogeneous 
sprays are seen to be very similar. The presence of chemical 
reaction decreases the radial velocity distribution for the 
homogeneous jets. Chigier and Roett (12) demonstrated that the 
velocity distributions may be expressed in the form 



• -r2 

2 c ^22 
I m 


(3.54) 


The spreading coefficients for the sprays described in Figures 11 and 
12 are shown in Table 1. The density variation between the injected 
fluid and the ambient fluid range from 1/1000 for the steam-water 
system to 660 for the oil drop-air system. The value of the spreading 
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Table 1 


Comparison of Jet Spreading Coefficients 


T3rpe of Jet Cj^ 

Isothermal air-air jet, Chigier and Roett (12) .0713 

Isothermal oil drop-air spray, Hetsroni and Sokolov (5) .0645 

Combusting gaseous jet, Chigier and Chervinsky (21) .0373 

Combusting spray, Chigier and Roett (12) .0620 

Air-water jet, Tross (8) . 0802 
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coefficient remains relatively constant for large density 
variations. 

Hetsroni and Sokolov (5) suggested tfie following expression for 
the radial spread of fluid concentration of a spray 


iy/Ay^ 


= Exp 




(3.55) 


where c^ is the spreading coefficient of the fluid concentration from 
the data of Hetsroni and Sokolov (5) 


c = .05 (3.56) 

m 

The value of Ay /Ay is fixed by the value of this quantity at 

o 

the penetration length position. Equation (3.55) is used to 
calculate the radial position of the boundaries of the flame, the 
liquid fuel, and the condensed water droplets. 


3.5 Evaluation of Integral Constants 

In order to complete the general model of the turbulent jet 
combustion process, the integral constants, entrainment constant, and 
unmixedness factor must be defined. Avery and Faeth (16) defined 
three constant groups from experimental data as follows 


14(2/12)^^^ E 
E(2/l2)^/^ I3 
I 


= 35.5 


0.096 


1 

E(2/l2)^/2 = 0.075 


(3.57) 

(3.58) 


(3.59) 
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1/2 

The constant group £( 2 / 12 ) =0.16 may be determined from the 

data of Ricou and Spalding (20) . 

Values of the integral constants may be determined by curve 
fitting velocity profiles from constant density air filter systems. 
The plot of combined velocity profiles is shown in Figure 13. The 
integral constants were determined to be equal to 

= 0.00631 (3.60) 

I 2 = 0.00319 (3.61) 

Table 2 summarizes the values for the thermal constants. 

Table 2 


Summary of Constants for the Variable Density Model 




Constant 

Value 



0.00631 

^2 

0.00319 

^3 

0.00383 


0.00383 

E 

0.00638 

e 

0.29 

r 
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Figure 13 Correlation of Jet Velocity Spreading Cliaract eristics 



CHAPTER IV 


RESULTS AM) DISCUSSION 


4.1 Introduction 

The major objective of the present Investigation was to examine 
the effectiveness of a locally homogeneous two— phase model for 
predicting the evaporation and combustion characteristics of a gas- 
liquid turbulent jet. It was assumed that the fuel droplets were in 
thermal equilibrium with the surrounding gas and moved at the local 
gas velocity. The liquid fuel was taken to be finely atomized at the 
injector exit and the flow was assumed to be fully developed 
immediately downstream of the injector exit. In addition the 
turbulent spray and flame boundaries are characterized by an Irregular 
wavy appearance as shown in Figure 8. Estimates of the spray and. 
flame boundaries represent the time-averaged position of the 
boundaries . 

A variety of test conditions were examined in the present study. 
Table 3 lists the test conditions for the type of boundary examined. 
The various properties used for the test fuels are listed in Table 4. 

4.2 Nonreacting Spray 

Experimentally determined radial spray boundaries for 

nonreacting pentane and methanol jets are shown in Figures 14 and 15. 

The edge of the flow field is also indicated on these figures by 

plotting the position where Ay/ Ay = ,01, for both the compressible 

o 

and incompressible models. The incompressible radial spread model 
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Table 3 


Siommary of Experimental Conditions 



Test Conditions^ 

Methanol 

Pentane 

Iniector Characteristics at Atmospheric 

Pressure 


Sauter Mean Diameter (ym)^ 

33 

28 

Maximum Droplet Diameter 

107 

79 

Spray Velocity of Injector Exit (m/sec) 

61.1 

68.5 

Cold Elow Tests 



Pressure (MPa) 

3, 6, 9 

3, 6, 9 

Length, of Flow Examined l-i— j 

*•0 •* 

90 

90 

Flame Boundary Tests 



Pressures (MPa) 

O.l,*^ 3, 6, 9 

0.1,^ 3, 6, 9 

Length of Flow Examined 

“0 

. 60 

90 

Spray Boundary Tests (Combustion) 



Pressures (MPa) 

3, 6, 9 

3, 6, 9 

Length of Flow Examined ['^j 

0 

Entire 

Boundary 

Entire 

Boundary 


a) injector pressure drop of 2.67 MPa. 

b) Calculated according to Muegle (25) . 

c) Entire flame boundary measured. 
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Table 4 


Summary of Physical Properties 


Property 

Methanol 

n-Pentane 

Water 


786.7 

626.14 


h^ KJ/kg 

fa 

-7640 

-2397.4 

-15865.8 

hfg KJ/kg 

-6276 

-2029.2 

-13422.9 

Sd I' 

2.505 

. 2.330 

4.18 

C KJ/kg K 
Pg 

1,370 

1.713 

2.15 

P (bar) 
c 

79.9 

33.74 

— 

T (K) 
c 

512.58 

469.6 

— 

A+ 

7.97328 

6.85221 

6.6788 


1515.14 

1064.63 

573.480 


232.85 

232.0 

260.0 


*^T in C and P in mmHg in Equation (3.42) 
^Evaluated at 1 MPa and 298 K. 





Figure 14 Nonreacting Spray 



for Methanol 
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closely approximates the compressible radial spread model in the 
case of the nonreacting sprays. 

In all cases the observed boundary of the spray lies well inside 
the conventional jet boundaries; the liquid does not appear to diffuse 
to the edge of the flow, as expected for locally homogeneous flow. 

The degree of radial liquid spread clearly increases as the pressure 
increases and the flow better approximates locally homogeneous flow 
at high pressures. 

Two effects could be acting to produce the wider predicted 
boundaries observed in Figures 14 and 15. First of all, the outer 
edge of the spray is composed of fine droplets, which may not be 
detected with the optical system used in the current investigation. 
Therefore, the complete width of the spray may be wider than observed. 

The second effect involves large drops in the spray following 
trajectories and not being diffused in the radial direction as 
required by a locally homogeneous model. For a given particle size, 
turbulent diffusion of particles is enhanced with increasing gas 
density. The test results, showing that the spread of the spray 
is greater at higher pressures, suggest that this effect is present 
to some degree. 

In a reacting spray, the drop size decreases with distance 
from the injector which helps to improve the approximations of the 
locally homogeneous model; results of this type are discussed in 


the next section. 
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4.3 Reacting Spray Results 

Predictions and measurements for methane and pentane sprays 
burning in air are illustrated in Figures 16 and 17, respectively. 

The results are given at pressures of 3, 6, and 9 MPa. 

At each pressure the experimental spray length is longer for 
methanol than for pentane. Predicted spray lengths also follow this 
trend, however, the predicted lengths are only about half as long as 
the measurements. As in the case of the nonreactive sprays, the 
radial spread of the spray Increases as the pressure Increases. 

In contrast to the nonreactlve sprays, there is a substantial 
difference between the compressible and incompressible predictions 
of the radial spread of the spray. The error in spray length 
influences the prediction of radial boundaries, and it is difficult 
to decide which model is best. In general, however, it appears that 
the compressible model overestimates the radial spread of the 
spray . 

The locally homogeneous models predicted that there was a region 
of the spray where water produced by the combustion process would 
condense. The predicted boundaries of this region, using the 
incompressible model, are also illustrated in Figures 16 and 17. In 
all cases, the liquid water boundary was contained within the liquid 
fuel boundary and the presence of water drops should not influence 
the experimental determination of the fuel spray boundary. 

For pentane, the model indicated that drop gasification was 
finally completed by the fuel passing through its thermodynamic 
critical point at pressures of 6 and 9 MPa. ,A11 other test conditions 
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involved conventional drop evaporation at a relatively constant wet- 
bulb condition. Comparison between theory and experiment in Figure 
17 shows no unusual change in the results for supercritical 
vaporization conditions. As in earlier studies of supercritical 
droplet combustion (1-4) , the total pressure must be greater than 
the critical pressure for supercritical gasification to occur; 
although specific limits were not determined for the present test 
conditions. The critical pressure of methanol (7.99 MPa) is too 
high for supercritical evaporation to be observed during the present 
tests. 

4.4 Flame Boundary Results 

Predicted and measured flame boundaries for methanol and 
pentane burning in air are illustrated in Figures 18, 19 and 20. 

The results are given at atmospheric pressure (0.1 MPa) and pressures 
of 3, 6, and 9 MPa. The entire flame boundary could be obtained at 
atmospheric pressure. At elevated pressures, the flames were too 
long to obtain the entire boundary with the present apparatus, and 
only the initial portion of the flow is illustrated. 

Flame boundaries at atmospheric pressure are illustrated in 
Figure 18. The locally homogeneous model gives a reasonably good 
prediction of the flame length for pentane, only slightly 
underestimating the measured length. Results are poorer for methanol, 
the predicted length is only about half as long as the measured value. 
Both models overestimate the radial position of the flame, however, 
the error is less for the incompressible model. 




Figure IS Flame Boundaries at Atmospheric Conditions 
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Figure 19 Flame Boundaries for Methanol 
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Results at elevated pressures are illustrated in Figures 19 and 
20. Similar to the spray boundaries, the measured width of the flame 
boundary increases as the pressure increases. The .compressible model 
continues to overestimate the radial spread, although the 
Incompressible model provides a reasonably good estimation of the 
predictions. The adequacy of the incompressible model must be treated 
with some reservation, however, since only the initial portions of 
the flame were observed at elevated pressure. 

4.5 Penetration Length Results 

The predictions of penetration lengths of sprays and spray flames 
were based on a correlation developed for gas-gas and gas-liquid 
systems (16) . The present measurements for liquid-gas systems are 
compared with the original data base, and the correlation, from 
Avery and Faeth (16) in Figure 21. The tolerances shown on the data 
represent the standard deviation for the test sample at each 
condition. A summary of the comparison between theory and experiment 
is given in Table 5. 

In general, the present measurements are above the correlation 
for both spray and flame boundaries. The model provides a qualitative 
estimation of the length of various phenomena, e.g., the flame length 
is correctly predicted to be substantially longer than the spray 
length, however, the theory underestimates the lengths by 30 to 50 
percent . 

Increasing the pressure causes a data point to shift toward the 
left-hand side of Figure 21. The spray length results indicate that 
measurements at high pressures are generally in poorer agreement with 




56 


10,000 r 


1,000 h 


100 t 


• 

i 

A 

Q 

O 

A 

O 

• 

o 

o 


PENTANE FLAME BOUNDARY 1 
METHANOL FLAME BOUNDARY i 
PENTANE SPRAY BOUNDARY T 
METHANOL SPRAY BOUNDARY’ J 
COj-Nj- SPRAY BOUNDARY, NEWMAN a BRZUSTOWSKI 
CI 2 -No -FLAME BOUNDARY, AVERY a FAETH 
35% Hj,* 65% CHjCHiCHjI 


CC>2 ♦ CITY GAS 

»Z 

CITY GAS 
CH CH CH 
CH'CH 
CO 

CITY GAS WOHL.ETAL 


FLAME BOUNDARY 
I HAWTHORNE, 

> ET AL 



CORRELATION, AVERY AND FAETH 
1/2 

« 35.5(m,/d„[T/>„M„) B) 


•KERNEY, ET AL 


I I I I I I t Ll I t I I I 1 I I 1 1 I L-I... L L l-U 1 1 1 — I I 1. 1 l-l 

1.0 10 100 1,000 

1/2 


Figure 21 Correlation of Penetration Length Data 
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Table 5 


Experimental and Predicted Penetration Lengths 


a 


Pressure 

(MPa) 

Type of 
Boundary 

Theoretical 

Length 

L/d 

o 

Experimental 

Length 

L/d 

o 

Percent 

Deviation 

Methanol 

0.1 

Flame 

3290 

6640 

51 

3 

Spray 

98.9 

142 

31 

6 

Spray 

70.9 

132 

47 

9 

Spray 

58.5 

104 

44 

Pentane 

0.1 

Flame 

6490 

7850 

17 

3 

Spray 

73.8 

132 

44 

6 

Spray 

52.8 

119 

56 

9 

Spray 

43.1 

100 

57 


^For combustion of sprays in air at 23 C. 





the predictions. Behavior of this tjrpe was not expected, since the 
locally homogeneous model should be a better approximation of the 
real flow as the density ratio of the two phases approaches unity. 

The spray length measurements of Newman and Brzustowski (10) 
for the system are also shown in Figure 21. Their measurements 

fall below the present correlation although they found good 
agreement with their model. It is difficult to explain this behavior 
in view of the present results. The measurements of Reference (10) 
were made close to the critical point and properties are difficult 
to estimate in this region, which provides one explanation. The 
fact that Newman and Brzustowski did not have to allow for 
unmixedness of the turbulent flow in order to achieve good agreement 
between their theory and experiment also suggests that this data might 
be atypical. No other result illustrated jLn Figure 21 could be 
correlated in this manner. Furthermore it is unusual for the 
penetration length of a liquid-gas system to be overestimated by the 
no-slip penetration length correlation. Any slip that occurred between 
the droplets and the surrounding gas would have the effect of 
lengthening the jet. 

4. 6 Distribution of Centerline Velocity 

The locally homogeneous Integral model. has been used to provide 
a correlation of centerline velocity using data for gas-gas and 
gas-liquid systems (16) , The correlation has the following form 

Mq . 1/2 

= .075 (TTp^M^) z/m^ 


moUc 


(4.1) 
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The comparison between Equation (4.1) and the steam-water jets 
measurements of Kerney, et al., (7), and the homogeneous jet 
measurements of Albertson, et al., (22) and Forstall and Gaylord (23) 
is illustrated in Figure 22. 

The correlation can be examined for liquid-gas jets using 
measurements provided by Newman and Brzustowski (10) for liquid 
carbon dioxide injected into mixtures of gaseous carbon dioxide and 
nitrogen at high pressures. The comparison is illustrated in Figure 
22. It is evident that good agreement is obtained between the 
correlation developed for gas-gas and gas-liquid systems, and the 
measurements of Reference (10) for a spray. 

4. 7 Discussion of Results 

Aside from the measurements of Newman and Brzustowski (10) , 
which seem to be atypical, the model consistently underestimates the 
length of spray and flame boundaries Two major factors could be 
responsible for this behavior: (1) inadequacies in the basic 

turbulent jet model, and (2) failure of the locally homogeneous model 
for the present two-phase flow. 

The advantage of the present turbulent jet model is that it 
provides simple analytical expressions for mean quantities in a 
variable density jet. While the model has been developed using a 
large data base for jets, it is not very sophisticated by present day 
standards. In particular, the development region near the injector 
is not treated explicitly, and the model is only adequate at large z/d 

0 

values. The effect of turbulent unmixedness, when treating a quantity 
such as a spray or flame length, is handled using a fixed value of e; 



M,/ ( rh„u 
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it is unlikely that a fixed value of £ is adequate for all 
circumstances . 

The potential effect of the oversimplifications of the 
turbulent model could be evaluated’ using a more sophisticated 
calculation. Khalil and Whitelaw (15) have recently reported results 
of this type. Their experiment involved combustion of a kerosene 
spray produced by a swirl atomized injector. Two injector conditions 
were examined, having cold flow Sauter mean diameters of 45 and 100 
ym, respectively. The injector was located at the axis of a swirling 
air flow. A locally homogeneous two-phase model was employed to 
analyze measurements of mean velocity and temperatures within the 
flame. A k-e turbulence model was employed in the flow calculations 
with a clipped Gaussian probability density distribution used to 
represent unmixedness in the reaction rate expression. 

Unfortunately, the base flow of Reference 15 is rather complex, 
which complicates the interpretation of the results. It was found, 
however, that the locally homogeneous model overestimated the rate 
of development of the flame, similar to the present findings. 

Results were poorer for the spray having the larger Sauter mean 
diameter. This suggests a progressive failure of the locally 
homogeneous model as the drop size increases. Both the present 
study and Reference (15) indicate that sprays having Sauter mean 
diameters as small as 30-45 ym are still only marginally represented 
by a locally homogeneous assumption. Based on this finding, it 
appears that weaknesses in the present turbulent flow model are not 
the major source of the discrepancies between measurements and 
predictions. 
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There are two aspects of the failure of a locally homogeneous 
model of a two-phase flow; (1) kinematic failure due to slip, and 
(2) failure of the assumption. of thermodynamic equilibrium. Slip 
limits the rate at which the jet momentum is transferred to the 
gas phase, reducing the entrainment of the ambient gas from the levels 
obtained for gaseous injection with the same total momentum in the 
flow. Therefore, a greater length is required to entrain 
sufficient ambient material to gasify or burn the fuel, as observed 
in the present experiments. 

Slip is particularly important for large drops, which carry a 
significant percentage of the momentum of the flow. Large drops also 
do not diffuse in a turbulent flow, which limits the radial spread 
of mass within a two-phase jet (26). Reduced levels of radial spread 
were observed in the present experiments for both combusting and 
noncombusting flows. 

The radial diffusion of particles approaches turbulent gas 
diffusion rates as the density of the two phases approaches unity 
(26) . The present measurements exhibit this tendency, with the high 
pressure sprays showing a greater lateral spread. However, penetration 
length predictions did not show a corresponding improvement. This 
behavior is due to the fact that penetration lengths of high pressure 

sprays are shorter, e.g., the theory implies that L is proportional 

- 1/2 

to p^ for a given injector condition. Therefore, the bulk of the 
process moves nearer to the injector at high pressures, where problems 
with slip are still important. 

The effect of the thermodynamic equilibrium assumption was not 
examined quantitatively during the present investigation, however. 
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it is clear that loss of thermal equilibrium due to the thermal 
inertia of the drops would tend to increase the length of spray 
and flame boundaries. 

In order to improve the model, a complete two-phase turbulent 
flow analysis must be considered. This Involves computing the life 
histories of the drops produced by the injector. The interaction 
of the drops and the gas phase must be considered, so that the 
environment of the drops can be determined. Analysis of this type 
have been reported by Crowe, et al., (27), and Jurewicz and Stock (28), 
although the results have not been confirmed by experiments. These 
calculations employ a k-e turbulence model to represent the gas 
phase, and allow for the distributed exchange of mass, momentum, and 
energy between the two phases. The most recent version allows for 
particle diffusion as well, although little data is available on 
particle diffusion rates for use in the calculations (28) . 

The difficulty in employing a comprehensive two-phase turbulence 
model for a spray is that a great deal of empirical information is 
required as input for the calculations. The drop size 
distribution of the injector, as well as the initial streamwise and 
radial velocity distributions for each drop size group, must be 
provided; this information is rarely available for practical Injectors 
under the hot firing conditions. There are significant uncertainties 
in other quantities used in the model: drop transport processes, 

evaporation, drag, etc., in a spray environment; turbulence modeling 
constants; turbulence generation and dissipation by drops; turbulent 

I 

drop diffusion; and proper representation of turbulent combustion 
processes (29) . 
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In view of all the difficulties- of a comprehensive model, the 
simplicity of the locally homogeneous flow assumption is very 
appealing. In view of potential uncertainties In input data and 
correlations used in a complete model, the accuracy of even the 
present simplified model may prove to be reasonably competitive 
(although the availability of a large number of parameters within a 
complete model provides greater scope for matching theory and 
experiment). Integral models, however, cannot be readily applied 
to the recirculating floxirs frequently encountered in practical 
combustion systems, and locally homogeneous models similar to the 
one used in Reference (15) are preferable in this case. 

The empirical constants used in the present model were fixed 
by earlier gas-gas and gas-liquid measurements (16) . Using these 
constants, the present test conditions left little room to adjust 
parameters and improve the predictions.. It was felt that the test 
range of this investigation was too limited to obtain adjusted 
parameters for sprays. The measurements employed a narrow range of 
injector velocities, only two fuels, and a single injector diameter. 
Additional data over a broader range of variables is necessary prior 
to attempting a generalized correlation valid for spray. 



CHAPTER V 


SUMMARY AND CONCLUSIONS 


5.1 Sijmmary 

This investigation considered a turbulent spray combustion 
process in a stagnant environment. Emphasis was placed on examining 
the capability of a locally-homogeneous two-phase flow model to 
predict the characteristics of the spray 1 The specific objectives 
of the study were as follows: 

1) Measure spray and flame boundaries for combusting and 

noncombusting pressure atomized sprays at various ambient 
pressures. Methanol and n-pentane were considered at 
pressures in the range 0.1-9 MPa. The higher pressures 
in this range exceed the thermodynamic critical pressure 
of the fuels. 

2) Available measurements were to be compared with predictions 
using a locally homogeneous flow model (16) . The model is 
based on measurements in gas— gas and gas-liquid jets, but 
had not been examined for sprays. 

The test configuration consisted of a single-hole, orifice-type 
injector, issuing into a stagnant environment of pure air. No swirl 
was applied to the liquid flow. The sprays had a Sauter mean diameter 
of 28-33 pm, under cold flow conditions, at atmospheric pressure. 

The major results of the investigation may be summarized as 


follows : 
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1) The length of the spray boundaries decreased and the radial 
spread increased with increasing ambient pressure. 

2) At comparable conditions, the spray length is shorter for 
pentane than for methanol. Flame lengths, however, show 
the opposite trend and methanol has the longer flame. 

3) The theoretical model, based on the locally homogeneous 
flow- approximation, correctly predicts the trends cited in 
Items 1) and 2) . The specific correlations are given by 
Equations (3.29) and (3.32). The correlations consistently 
underestimate the measured lengths, with errors on the 
order of 30-50 percent. Although the present turbulence 
model is relatively crude, similar behavior has been 
reported for a k-e turbulence model of a spray, using the 
locally homogeneous flow approximation (15) . 

4) ,The best prediction of radial boundaries was obtained using 

radial length scales for incompressible flow. This agrees 
with other measurements of variable density flows where it 
has been found that spread rates are similar to constant 
density flows (5, 13). 

5) Mean velocities in a spray were predicted reasonably well 
with the present model for z/d^ > 10. 

6) The accuracy of the predictions is poorer at high pressures 
even though the density ratio of the two phases approaches 
unity which should tend to improve the locally homogeneous 
approximation. This behavior appears to be due to the fact 
that the penetration length decreases as the pressure 
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Increases, therefore, the bulfe. of the process approaches 
the injector where slip effects are still important. 

7) The predictions indicate the presence of regions within a 
burning spray where water vapor produced in the combustion 
process should condense. For the conditions of the present 
experiments, the liquid water boundary always fell within 
the spray boundaries. 

8) Predictions for pentane at 6 and 9 MPa indicate that the 
spray finally gasified by the drops passing through the 
thermodynamic critical point. The results do not indicate 
any unusual phenomena when this occurred. 

5.2 Conclusions 

The conclusions of the investigation are as follows: 

1) The present locally homogeneous spray model is convenient 
to use, and requires a minim um amount of input data, but 
tends to overestimate the rate of development of the 
combustion process. Quantitative predictions could be 
improved by adjusting some of the empirical parameters 
within the model, from the values that were optimized for 
gas-gas and gas-liquid jets. However, the present test 
range was felt to be too limited to undertake an adjustment 
of this type, A larger data base, particularly including 
more fuel types, different injector velocities and 
different injector diameters is necessary before this step 


can be undertaken with some confidence. 
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2) The simple integral model uaed to represent th.e turbulent 
jet in this- investigation is seriously limited in its 
capabilities to be extended to recirculating flows. Further 
work should consider more complete turbulence models, e.g., 
the model used byRhaliland Whitelaw (15), which are capable 
of this extension, A model of this type would remove some 
of the uncertainties in evaluating the locally homogeneous 
flow assumption. Efforts to date in two-phase flows have 
been limited, and further study is required to establish 
the range of validity of the locally homogeneous flow 
approximation for spray evaporation and combustion processes. 

3) Spray combustion predictions attempted thus far have 
employed the locally homogeneous flow approximation. Models 
of this type have consistently overestimated the rate of 
development of the combustion process. Mor.e exact results 
require the development of complete two-phase flow models. 
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APPENDIX 


EXPERIMENTAL RESULTS 


Table 6 


Spray Boundary Data, Noncombusting Methanol Jet 



Pressure : 

3 MPa 

6 MPa 

9 MPa 

z/d^ 

t/To 

2/‘io ^/^O 

z/do r/r^ 


10.3 

1.9 

9.7 

2.9 

10.0 

3.8 

20,6 

3.7 

20.3 

5.0 

20.4 

6.7 

30.9 

4.5 

30.1 

6.7 

30.2 

9.5 

41.2' 

6.7 

46.5 

8.6 

40.2 

11.4 

51.6 

8.3 

50.8 

11.3 

51.3 

16.2 

61.9 

11.1 

61.2 

14.9 

61.3 

19.9 

71.5 

13.3 

71.4 

16.8 

71.4 

22.3 

82.5 

18.0 

82.3 

19.8 

81.8 

25.0 
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Table 7 


Spray Boundary Data, NoncomBustlng Pentane- Jet 



Pressure: 

3 MPa 

6 MPa 

9 MPa 

z/do 

rfr 

o 

z/d^ r/r^ 



10.1 

3.9 

10.2 

4.9 

10.0 

3.5 

20.6 

6.2 

20.1 

5.6 

20.0 

6.6 

29.9 

7.9 

36.4 

8.4 

30.0 

8.0 

41.6 

9.6 

41.0 

12.0 

40.2 

12.9 

51.8 

13.0 

51.2 

14.9 

51.5 

15.2 

61.4 

14.9 

61.1 

16.7 ■ 

61.6 

17.6 

71.6 

15.6 

71.4’ 

19.0 

71.4 

19.9 

82.3 

17.5 

81.6 

20.2 

81.8 

21.1 

92.7 

18.5 

91.8 

21.3 

92.8 

• 23.5 
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TaBle 8 


Spray Boundary Data, Combusting Metlianol Jet 


Pressure 

5 .. 3 MPa 

6 MPa 

9 MPa 

z/do 

-/-o 

z/do 


z/d^ 

r/ro 

11.4 

3.8 

13.3 

5.0 

13.1 

12.3 

43.0 

8.2 

43.0 

10.2 

42.8 

11.5 

71.8 

12.5 

73.1 

13.6 

72.9 

14.6 

104.9 

12.6 

104.0 

15.1 

104.2 

0.0 

142.3 

0.0 . 

132.4 

0.0 ■ 




Table 9 


Spray Boundary Data, Combusting Pentane Jet 


Pressure; 

3 MPa 


6 MPa 


9 MPa 

z/d 

o 

r/r 

o 

N , 

D. 

O 

r/ro 

z/do 


30.0 

6.7 

30.0 

9.3 

30.0 

9.4 

60.0 

15.2 

60.0 

14.0 

60.0 

20.1 

90.0 

16.1 

90.0 

20.2 

90.0 

23.6 

112.9 

15.5 

119.3 

0.0 

10.0 

0.0 

131.8 

0.0 
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Table 10 


Flame Boundary Data, Methanol Jet 


Pressure ; 

0.1 MPa 

3 MPa 

6 MPa 

9 MPa 


z/d 

o 

r/r 

o 

z/d 

o 

r/r 

o 

z/d 

o 

r/r 

o • 

z/d 

o 

r/r 

o 

14.9 X 10^ 

115 

15.0 

12.8 

15.0 

14.0 

15.0 

15.6 

30.0 X 10^ 

135 

45.0 

22.2 

45.0 

37.3 

45.0 

29. -2 

49.8 X 10^ 

203 

60.0 

25.5 

60.0 

32.8 

60.0 

34.9 

59.8 X 10^ 

155 







6.64 X 10^ 

0.0 
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Table 11 



Plame Boundary- Data, 

Pentane 

Jet 



Pressure: 0.1 MPa 

3 MPa 

6 MPa 

9 MPa 

o 

'll 

N 

r/r^ 

z/do 

r/r^ 

z/do- 

r/r 
' o 

r/do 

r/r^ 

.89 X 10^ 

246 

30 

10.4 

30 

13.4 

30 

14.5 

1.59 X 10^ 

433 

60 

18.5 

60 

24.7 

60 

25.8 

1.63 X 10^ 

588 

90 

25.3 

90 

40.6 

90 

41.7 


3.44 X 1Q3 650 

4.34 X 10^ 763 

5.35 X 103 700 

6.18 X 10^ 525 

7.06 X 10^ 338 

7.85 X 1q3 


0.0 
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